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ABSTRACT 

Generalized three body problem is investigated with one of the point masses replaced with 
spherical body of finite size. The model provides simplified description of the Thorne-Zytkov 
objects, protoplanets, red supergiants and other core-halo structures in binary systems. 

Proposed description closely resemble the planar restricted circular three body problem. 
Mechanical system includes two point masses and third spherical body with prescribed den- 
sity. Only purely gravitational interactions are considered. System is shown to be unstable due 
to parametric resonance between orbital and internal frequencies. The central compact object 
can be ejected if resonance occur during the evolution of the system. Some interesting con- 
clusions from this discovery are presented. In particular, these findings may be important for 
core-collapse supernova theory, because helium core might be placed off-center relative to the 
hydrogen envelope. The instability is also possible during protoplanet growth, with hypothet- 
ical ejection of the rocky core from gas giants and formation of the "puffy planets". For sake 
of curiosity I note, that the Earth-Moon, and the Earth-Sun systems are stable in the above 
sense, so no hope for mankind exist in case of artificial black hole formation in the course of 
e.g. LHC experiment. 

Key words: Physical data and processes: gravitation, instabilities, chaos — Astrometry and 
celestial mechanics: celestial mechanics — Planetary systems: planets and satellites: dynam- 
ical evolution and stability, formation, planet-star interactions — Stars: binaries: general, su- 
pergiants, kinematics and dynamics 



1 INTRODUCTION 

Long established approach to astrophysical binaries and related dy- 
namics assume that mechanical system can be described in terms 
of point masses. However, in many situations, this assumption is 
clearly broken, with Thorne -Zytkov |Thorne & Zytkow|1977) ob- 
ject being extreme e xample l|Thorne & Zytkow|1977 Vanture et al.| 



p99l|Cannon|1993||Cannon et al.|1992J . Thorne-Zytkov object is 
composed of ordinary star and compact object (either black hole or 
neutron star) in the very center. Sometimes one consider other simi- 
lar situations e.g. non-standard Solar models with black hole inside 
Piawking|197T] |Clayton et al.|1975) , or artificial one created on 
Earth in failed high-energy experiment ( [Casadio et al.|20lO{ |Ble-| 
icher & Nicolini 2010, Gingrich 20101. To less extent, usual red 
giant pre-supernova stars are also more point-mass core plus ex- 
tended halo than single body, see e.g. Wo osley et al.| ( |2"002} . An- 
other example where presented model might be useful is dynam- 
ics of icy moons (e.g. Europa, Enceladus) and exoplanets („blue 
ocean" super-earths Haghighipour (20111), with rock or ice core 
surrounded by extended fluid ocean ( Kuchner 2003 , Leger et al. 
|2004l >. 
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Main motivation for this model is suggestion of |Arnett &| 
Meakin (2011), that in core-collapse supernovae iron core might 
be displaced with respect to geometrical symmetry center of the 
extended and usually more massive hydrogen envelope. Arnctt & 
|MeakinH2011) proposed hydrodynamical L=l instabilities during 
Si burning as a main cause of the displacement. Here I propose 
another model: gravitational instability in binary system^] 

Goal of this work is to answer the following question: is the 
point mass in the center of extended object stable with respect to 
perturbations caused by the third body orbiting the pair? Naively 
one might quickly answer „no" because it seems always possible to 
tune internal frequency (e.g. changing equation of state) and orbital 
frequency (e.g. by inspiral or adding mass to the system). How- 
ever, many years of struggle to answer very similar questions in the 
classic three body problem ( Mardling 2008a), as well as many sim- 
plifying assumptions, suggest caution and rigorous mathematical 
approach. 



1 |Arnett &~M eakin ( 2011) wrote: ,, If there were a driving mechanism 
for core-mantle oscillation, here would be an asymmetry due to the dis- 
placement of the core and mantle relative to the center of mass." Presented 
dynamical model describes such a driving mechanism in close binaries. 
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2 RESTRICTED PLANAR CIRCULAR THREE BODY 
PROBLEM APPROACH 

Mechanical system of interest includes two masses m and M or- 
biting the center of mass on circular orbits. Mass m (first body) is 
a point mass. Mass M (second body) is extended spherical body 
with density p(r). Rotational degrees of freedom for mass M are 
not considered, and body is rigid. Third body is a test body, and 
its mass is assumed to be negligible (see Sect. [4] for more general 
model). In this section I further assume that density p(r) inside 
mass M is constant to simplify formulae. Additionally, mass m is 
assumed to orbit outside radius R of mass M . Distance between 
geometrical center of mass M and m is equal to d, i.e., d ^ R. 
Third test body is initially in the center of mass M, where also co- 
ordinate origin is placed. In co-rotating system, equations of motion 
for third body, restricted to orbital plane, are: 



x — 2 u) y + k x+ 
y — 2 cj x + k y+ 
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System ([TJ is very similar to the classical planar restricted cir- 
cular three body problem, see Capiski & Zgliczyski (2011). Note, 
that system {T}, without terms explicite involving G, describe Fou- 
cault pendulum problem ( Landau & Lif shitzJT969j . 

Conserved „energy" for system |TJ is: 



E=-x 



' ' 2 + \v +U(x,y), 



System is considered unstable with respect to small perturba- 
tions if at least one solution of (7} has a positive real part: 

Re(X) > 0. 

Resolving above conditions lead to the instability criteria: 
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Criteria ([8} has been verified solving {T} numerically, and 
solving linearized system {5} analytically. Result {8a| can be ob- 
tained from analysis of the potential \4b\ extremum at x = 0, y — 
as well. Left-hand side of ^8aJ, i.e., condition 4/37rp > M/d 3 , 
from astrophysical point of view it is not interesting because the 
central density cannot be smaller than mass divided by volume 
of the ball. This if especially true for any spherical astrophysical 
body of interest, because density always decrease outwards. The 
same argument apply to ( |8b[ l. Even in the most favorable situation 
m — 2M, the central density should be less than half of the average 
density for mass m for this instability to occur. 

Therefore we may simply write simplified form of d8p, rele- 
vant to the astrophysical applications: 
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Note that system {TJ and energy \Ab\ are suitable for small pertur- 
bations of the test body only. To explore possible ejection of the 
test body from mass M we have to explicitly handle p(r), or at 
least consider radius of mass M. This is done in Sect. [3] 

Linearization of the equations {T} for small perturbations 
around point x = 0, y — lead^to the linear system: 



C-2aj£+(fc-2g)C = 0, 



^ + 2cu( + {k + q)£ = 0, 



where: 
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Eigenvalues A of the system |5j are solutions of the algebraic equa- 
tion: 



(A 2 + k - 2<7)(A 2 + k + q) + Alu 2 = 0. 



(7) 



3 NUMERICAL VERIFICATION AND THE 
LONG-TERM BEHAVIOR 

System {T} is relevant only to stable situations, and is a good tool 
to derive the instability criteria (8j. For unstable cases point mass 
is likely to abandon central region. Actual density distribution p(r) 
(location of „surface radus" defined as p(R) — in particular) is 
essential. We have to replace harmonic oscillator (,, spring") force 
with more realistic one. Introducing function /: 



m = -g 



the system of equations of motion is similar to l|T|(, but now: 

k = fi^x 2 + y 2 ) - uj 2 . (10) 
Mass M is now equal to mass inside distance d: 



2 After substitution x(t) = e £(t),y(t) = e £(i) into {TJ, series expansion 
has been calculated with respect to e, and higher-order terms dropped. 
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and we must further assum^] that „radius" R < d, that is, area 
where the density distribution p(r) > is encircled by distance d. 
We still are able to provide conserved energy: 
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where r = y/x 2 + y 2 and: 
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Constant C depends on lower integration limit R, and should be 
used to remove potential singularity at r — 0. 

Global dynamics is easy to predict qualitatively using en- 
ergy considerations. Three typical cases are presented for uniform 
sphere of density p and radius R (dashed circle) in Fig.[TJ Allowed 
„Hill region" has been shaded. In stable situation (Fig. [TJ panel 
(a) ), we have three disconnected regions: central area of mass M, 
neighbourhood of mass m, and ,, outer space" extending to the in- 
finity. Perturbed test body simply oscillates with frequency related 
to the „central" density, simultaneously rotating like Foucault pen- 
dulum. In presence of drag it will settle down at the geometrical 
center of mass M if perturbed. When the density drops to value 
below critical l[8j, central region of mass M and area surround- 
ing mass m become connected, cf. Fig.[TJ middle panel (b). Test 



3 If mass m is inside region where p(r) > the inertia! mass Mi nert is 
different than the gravitational mass M grav , and equal to: 

poo 

M inert = 4tt / p(C)C 2 ^C 
Jo 

while the gravitational mass is: 
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Therefore, for mass m inside density distribution, circular two-body prob- 
lem has slightly different solution compared to ,, classic" one. Kepler fre- 
quency is: 

Mgrav + m Mi™"t 

Unfortunately, third test body still have equal gravitational and inertial 
mass, and therefore do not co-rotate with mass M\ 



G- 



body begin to oscillate with growing amplitude forced by gravita- 
tional pull of mass m. After some time, depending on amplitude 
of initial perturbation, test body is ejected from mass M and en- 
ters chaotic orbit around mass m, see Fig. ^ panel (b). It is still 
bounded with binary system. With density much smaller than criti- 
cal, all three allowed regions become connected, and the test body 
will be ejected from the system, spiraling out into infinity (Fig.[T] 
panel (c)). Described behavior has been confirmed by the numerics. 
Example trajectories are shown in Fig.[T]using red curves. 

Analyzed instability is not strictly resonant behaviour with 
Kepler frequency {2), because instability criterion |8a| includes 
additional factor 3 before mass m. However, using language of 
resonance we can easily understand underlying physics (Mardling 
2008b). Consider classical textbook example of uniform density 
ball with radius R and satellite in circular orbit just above surface 
(d = R). Well known result is that internal harmonic oscillator and 
orbital frequencies are identical. Therefore, test mass in the center 
is in resonance with small mass in orbit. However, astrophysical 
bodies are centrally condensed (p c > p) and binary companions 
are usually more distant (d > R), so this resonance disappear. But 
orbital frequency depends on both masses M and m. Adding mass 
to the system in form of larger mass m we can make these fre- 
quencies equal again. In practice, for stellar binaries, variables are 
central density and/or orbital distance, while the total mass of bi- 
nary is conserved. Note, that we have never assumed that mass M 
is dominant. In fact, instability is more likely if mass m ^> M. 

Inevitable presence of ,,drag" and „friction" forces acting on 
a test body inside mass M has only a minor effect on result- 
ing dynamics. It may suppress amplitude of oscillations and de- 
lay onset of the instability, but due to ,, top-hill" position of point 
x = 0, y — in l |4b[ ) below critical central density, the global dy- 
namics remains the same. Numerical tests with drag of the form: 



K(p)r 



(13) 



revealed, that main effect is to delay onset of the instability un- 
til ejection from the system becomes possible, cf. Fig. [T] panel 
(c). Therefore, ejection is more likely outcome, while chaotic or- 
bits (Fig.[TJ panel (b)) might be rare. After test body leave mass M 
evolution is identical to the classical restricted planar circular three- 
body problem. With p(r) — > S(r) mass M becomes point-like, and 
we get back to the classical three body problem. 
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4 THREE BODY INSTABILITY 

In this section we remove assumptions of planar motion and neg- 
ligible mass of the third body. Full system of nine equations de- 
scribing motion of the three masses without additional simplify- 
ing assumptions can be transformed to the co-rotating system as 
well. Very lengthy calculations show that criteria {8} nearly survive. 
Solving linearized eigensystem we are able to show, that instability 
is present if: 



M 4 



M + 3m(l + f i/At)- 1 

J 3 



(14a) 



In the above example „unstable core" has a minimal mass of ~ 
4.3 Mq. Noteworthy, edge of the He core is placed at ~ 4.2 Mq. 
This result do not significantly vary within mass range 1 Mq < 
m < 100 A/0. I conclude that splitting the red giant into helium 
core and hydrogen envelope is the most appropriate. If we treat He 
core as a point mass, model presented in the article can be applied. 
It is likely, that during supernova event in close binary system, ex- 
plosion engine will be displaced with respect to geometrical cen- 
ter of the hydrogen envelope. Recently discovered stripped helium 
core (low-mass white dwarf) in binary system (Maxted et al .|201 1| | 
might have been formed similar way. 



4 lmM + ri — m/8 .„ „, s 
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New factor in \14\ is the mass of the third body equal to fj,. 

Again, only \l4a\ is of astrophysical interest. Using orbital 
frequency: 



2 G(m + fi + M) 

uj — — 



and internal frequency: 



4 n 



we may write j!4a| > as: 

2 G{m + n) 2 2. G(m + (i) 3Gm /A ^-i 
w + Js <uj c <u+ ^ + _(l + M /M) . 

Instability is a consequence of internal and orbital frequency over- 
lap. Width of the resonance is proportional to „forcing" mass m, 
and reduced by factor dependent on mass ratio (j,/M. The most 
important factor is mass m, because (1) increase orbital frequency 
allowing for resonance (2) increase width of the instability window. 



5.2 Exoplanets and moons 

Exoplanet formation and structure often consider compact core ac- 
creting mass in the form of extended low-density envelope (Laugh- 
|lin et al.|2004| >. In this situation, dominant is the mass m, i.e., the 
central star, and instability occur if: 



T 



1 g/cc \1 day 



< 0.057 



Alt, 



(17) 



where: p is the density of the envelope, M env is the mass of ac- 
creted envelope and Mtot is the total (core+envelope) mass of the 
protoplanet. 

It is not surprising that all of the tested exoplanets are sta- 
ble according to criterion (17) . This is also true for so-called Ultra 
Short Period Planets ( Sahu et al. 2006 1 with orbital period less than 
a day. However, stability margin is often small. We may specu- 
late that some of the „puffy planets" (Hart man & et. al.|2011) , i.e., 
very low density lupiter-like objects close to the central star, were 
formed in process involving ejection of the dense planetary core 
due to instability presented in Sect. [4] Even if the instability do not 
lead to the core ejection due to e.g. friction, dissipated energy might 
inflate the planet. This, however, does not explain lack of the rocky 



5 ASTROPHYSICAL EXAMPLES OF THE INSTABILITY 
5.1 Massive binaries 

Application of the results from Sect. [2] to a massive star is not 
straightforward, because ,,core" is not well-defined and separated 
from envelope. Red supergiants are indeed objects with nearly 
point-like core and extended low-density envelope. However, split- 
ting radius is more or less arbitrary. It is also not clear what re- 
ally will happen if instability is operational. Thorne-Zytkov objects 
are exception, because central object is well approximated by point 
mass. 

To overcome mentioned difficulties, I adopted the following 
procedure: star with total mass M» and radius R, is artificially 
divided into two parts: (i) central „core" region with r < £ (ii) 
outer envelope with R, > r > £. Now, instability criteria \14\ are 
functions of parameter £ with: 



H = fh(£), M = M„-m(£), 



(15) 



where rh(£) denotes mass enclosed by sphere with radius £, To 
further reduce complexity (dimensionality) of the analysis I as- 
sume that perturbing mass m is as close to star as possible, i.e, 
d = R„. Stellar model sl5 of |Woosley et al.| (2002) with mass 
M* = 12.8M and radius R* = 3.85AU at Si burning stage has 
been used as an example. The system is stable if: 



-TTp(0Rl > (M, + 3m) 1- 



m(0 



(16) 



5.3 LHC black hole 

It has been speculated that LHC or other high-energy experiment 
might produce artificial black hole, that do not explode immediately 
vis Hawking radiation (Casadio e t al. 2010; Bleicher & Nicolim] 
|2010| [Gingrich 20T0)l. Such a black hole would settle at central re- 
gion of the Earth and consume the planet. I have applied instability 
criterion |9]( to the Earth-Moon and Earth-Sun systems. Unfortu- 
nately, central position is stable by a wide margin in the sense of 
instability J8l. 



6 CONCLUSIONS 

Presented generalized three body model has been analyzed using 
analytical and numerical techniques. Instability of the point mass 
placed in the geometrical center of mass was found. Analytical cri- 
teria were derived using linearized system, and verified numeri- 
cally. Proposed description has been assumed to be valid for as- 
trophysical binaries, where one of the companions has a core-halo 
structure. At least two possible sites for the instability were found 
: massive red supergiants in binary system (Subsect. \5A\ and for- 
mation of the exoplanets (Subsect. |5. 2\ . Model can be applied also 
to a more exotic situations, e.g., Thorne-Zytkov objects or central 
black holes of astrophysical and artificial origin. 
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Three body model need to be validated. In particular, astro- Woosley S. E., Heger A., Weaver T. A., 2002, Reviews of Modern 

physical bodies of interest are not rigid (typically gaseous), and Physics, 74, 1015 

might not react to a driving force as a whole. Orbital timescale 
Torb ~ -R 3//2 is however much longer compared to sound crossing 
time Tsnd ~ R, at least for red supergiants with R 3> Rq : 



T -^0.15x-/ M/Mg 



'Torb y "7 

Tq R/Rq 

where 7 is the polytropic exponent, Tg is temperature divided by 
10 9 K. Ultimately, three dimensional hydrodynamic model with 
appropriate treatment of external gravity source, either analytic or 
numerical, should be used to verify instability. 

Binary and multiple systems in astrophysics are rather a rule 
than exception, as well as core-halo structure of components, in- 
cluding dark matter halos with central black hole as well. There- 
fore, instability presented in the article might be of a very common 
occurrence in the nature, influencing formation and evolution of the 
astrophysical bodies and structures on various scales. 
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